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Abstract 

O ■ 

Lie and Q-conditional symmetries of the classical three-component diffusive Lotka- 
W ■ Volterra system in the case of one space variable are studied. The group-classification 

t"""- problems for finding Lie symmetries and Q-conditional symmetries of the first type are 

completely solved. Notably, non-Lie symmetries (Q-conditional symmetry operators) for 
a multi-component non-linear reaction-diffusion system are constructed for the first time. 
An example of non-Lie symmetry reduction for solving a biologically motivated problem 
is presented. 

> \ 
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1 Introduction 

In 1952 A.C. Turing published the remarkable paper pQ, which a revolutionary idea about 
mechanism of morphogenesis (the development of structures in an organism during the life) 
has been proposed in. From the mathematical point of view Turing's idea immediately leads to 
construction of reaction-diffusion systems (not single equations !) exhibiting so called Turing 
instability (see, e.g., Chapter 14.3 in [2]). It should be stressed that nonlinear react ion- diffusion 
(RD) systems are govern equations for many well-known nonlinear second-order models used to 
^ describe various processes in physics [3 J , biology (2J H] and ecology [5] . Since 1952 nonlinear RD 

systems have been extensively studied by means of different mathematical methods, including 
group-theoretical methods. Nevertheless finding Lie symmetries of two-component RD systems 
was initiated about 30 years ago [6], complete solving this group-classification problem was 
finished only a few years ago in papers [7j |8] (for constant diffusivities), [9j [lOj [TT] (for non- 
constant diffusivities) and [12J (for constant cross-diffusion). It should be stressed that there 
are only a few papers devoted to search for non-Lie (conditional, non-classical) symmetries of 
such systems [131 HH US] • Because finding non-Lie symmetries for RD systems is very difficult 
problem [T4"| [16], only special cases of the general two-component RD system were examined 
in these papers. 

All the papers cited above, excepting [TT], deal with the two-component RD systems. In 
section 4 of [TT], Lie symmetries of a class of multi-component RD systems are completely 
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described. To the best of our knowledge there are no any papers devoted to construction of 
conditional symmetries for the multi-component RD systems. On the other hand, nonlinear 
mult i- component RD systems is an important tool for mathematical modeling a wide range 
of processes involving several kinds of species (cells, chemicals etc.). Moreover, such systems 
possess some properties, which are not common for relevant two-component systems. Thus, it 
is time to extend the results obtained for two-component RD systems on the multi-component 
systems. It turns out that this is a highly non-trivial problem and this paper is devoted to 
solving this problem for the well-known three-component RD system. 

In the present paper, we shall consider the diffusive Lotka-Volterra (DLV) system 

Ai tit = u xx + u(ai + biii + C\V + d]w), 

X 2 v t = v xx + v(a 2 + b 2 u + c 2 v + d 2 w), (1) 
X 3 w t = w xx + w(a 3 + b 3 u + c 3 v + d 3 w), 

where u(t, x), v(t, $), w(t, x) are unknown concentrations, a^, df. and A& > are arbitrary 
constants (hereafter k — 1, 2, 3 and the subscripts t and x denote differentiation with respect to 
these variable. These constants have the relevant biological (chemical) interpretation depending 
on a type of interaction between populations (cells, chemicals), which system (TjQ) describes 
[21 HI [5]. Setting formally w = 0, we obtain the two-component DLV system. Lie symmetries 
of such system have been completely described in [20] , while our recent paper [15] was devoted 
to search of its Q-conditional symmetries. 

System (pQ) is the standard generalization of the classical Lotka-Volterra system that takes 
into account the diffusion process for interacting species (see terms u xx , v xx and w xx ). Never- 
theless the classical Lotka-Volterra system was independently introduced by A. J. Lotka and V. 
Volterra about 90 years ago, its different generalizations are widely studied at the present time 
because of their importance for mathematical modeling many processes in physics, biology, ecol- 
ogy etc. Notably, there are several recent papers devoted to rigorous study of multi-component 
diffusion Lotka-Volterra system (see [T7 1 fT8 t fT9] and papers cited therein). 



2 Main Definitions 

Because the DLV system ([T|) belongs to the general class of three-component RD systems 
Ai«t = u xx + C 1 ^,^, w), 

X 2 v t = v xx + C 2 (u,v,w), (2) 
X 3 w t = w xx + C 3 (u, v, w), 

where C h (u,v,w) are arbitrary smooth functions, we formulate the main definitions for any 
system of the form (J2]). 

To find Lie invariance operators, one needs to consider system (T5]) as the manifold Ai = 
{S 1 = 0,S 2 = 0,S 3 = 0}, where 

Si= Xiu t - u xx - C l {u,v,w), 
S 2 = X 2 v t - v xx - C 2 (u,v, w), 
S 3 = X 3 w t - w xx - C 3 (u, v, w), 
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in the prolonged space of the variables: t, x, u, v, w, u t , v t , w t ,u x , v x , w x , u xx , v xx , w xx , u xt} v xt , w xt , 
u tt , Vtt, Wtt- According to the well-known definition, system (|2j) is invariant under the Lie group 
generated by the infinitesimal operator 

Q = x -> u i v i U ')^t+C 1 (^ x i u i v i w )d x +ri l (t, x, u, v, w)d u +r] 2 (t, x, u, v, w)d v +i] 3 (t, x, u, v, w)d w 

(3) 

if the following Lie's invariance conditions are satisfied: 





C l (u, v, w)) 




o, 




M 




V xx 


C 2 (u, v, w)) 




0, 




M 






- C 3 (u, v, w)) 


M 


= 



(4) 



Q(S 3 ) = Q(X 3 w t 



where the operator Q is the second prolongation of the operator Q (see, e.g., [2H [221 1231 121]). 

Definition 1. [16J Operator ([3]) is called the Q-conditional symmetry of the first type for the 
RD system ([2]) if the following invariance conditions are satisfied: 



Q(Si) = Q(\\u t - u xx - C l {u, v, w)) 

2 2 ' 

Q(S 2 ) = Q(X 2 v t - v xx - C 2 (u,v,w)) 

2 2 v ' 

Q{S 3 ) = Q(X 3 w t - w xx - C 3 (u, v, w)) 



Mi 



Mi 



o, 

= 0, 



(5) 



Mi 



where the manifold Aii is either {Si = 0, S 2 = 0,S 3 = 0,Q(u) = 0}, {Si — 0, S 2 — 0,S 3 = 
0,Q(v) = 0} or {Si = 0,S 2 = 0,S 3 = 0,Q(w) = 0}. Hereafter the notations Q(u) = 
£,°u t + ^u x — r] 1 , Q(v) = £°i>i + ^v x — rf and Q(w) = £,°w t + — rf are used. 
Definition 2.[16j Operator ([3]) is called the Q-conditional symmetry ( non-classical symmetry 
in terminology used, e.g., in [23]) for the RD system ([2]) if the following invariance conditions 
are satisfied: 



Q(Si) 

2 


= Q(XlU t - u xx - 
2 


- C l (u, v, w)j 


= 0, 

Ms 


Q(s 2 ) 

2 


EE Q(\ 2 V t - V xx - 


- C 2 (u, v, w)) 


= 0, 

Ms 


Q(s 3 ) 

2 


EE Q(X 3 w t - W xx 


— C 3 (u, v, w)) 


= 

M :i 



where the manifold M 3 = {Si = 0,S 2 = 0, S 3 = 0, Q(u) = 0, Q(v) = 0, Q(w) = 0}. 

It is easily seen that Ai 3 C Mi C M, hence, each Lie symmetry is automatically a Q- 
conditional symmetry of the first type, while each Q-conditional symmetry of the first type 
is also Q-conditional (non-classical) symmetry (see also [16] for extensive discussion about a 
hierarchy of conditional symmetry operators). 
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Proposition 1. Let us assume that X\ = h 2 (t, x,u,v,w)d v + h 3 (t, x,u,v,w)d w (hereafter 
h 2 and h 3 are the given functions) is a Lie symmetry operator of the RD system ([2]) while Qi is 
a Q- conditional symmetry of the first type, which was found using the manifold Ai\ = {Si = 
0, S 2 = 0, 5*3 = 0, Q{u) = 0}. Then any linear combination C\X\ + C 2 Q\ (hereafter C\ and 
C 2 are arbitrary constants) produces another Q-conditional symmetry of the first type. 

It should be stressed that this statement is not valid for arbitrary given Q-conditional 
symmetry but only for that of the first type. Proposition 1 will be applied in the next section 
to prove Theorem 2. 



Obviously, the RD system (CQ) for arbitrary functions C k (u,v,w) admits the two-dimensional 
Lie algebra, called the principal (or trivial) algebra, with the basic operators : 



It can be easily shown that ([6]) is the principal algebra also for the DLV system (CQ). Note that 
we want to exclude the semi-coupled systems, i.e. those containing an autonomous equation, 
hence, hereafter the restrictions 



are assumed. 

To find all possible extensions of principal algebra in the case of the DLV system (TjQ), one 
needs to apply the invariance criteria (HI) and to solve the obtained system of determining 
equations (DEs). Omitting rather standard calculations we present the system of DEs in 
question: 



3 Main Results 



p t = d t , p x = d x . 



(6) 



(7) 



f° = f° = f :0 = f = p 1 = f 1 = f 1 = n 

>M >>V ^>W "sMI U 5 



(8) 
(9) 
(10) 

(11) 
(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

(22) 



Vuu Vuv Vvv Vww Vuw Vvw "J ™ 1) ^) 




(Ai - X 2 )vl = (Ai - A 3 )^ = 0, 
(Ai - X 2 )vl = (A 2 - X 3 ) V 3 V = 0, 
(Ai - XsK = (A2 - A 3 )^ = 0, 
2£ - # = 0, 

2£%L + (A 2 - X^vl = 0, 

^°vL + (A 3 - Aoe 1 ^ = 0, 

Hu + Ai^ 1 = 0, 

+ A 2 et = 0, 

*vL + As^ 1 = 0, 

V'CI + r, 2 Cl + V 3 Cl + ri x - X lV l + (2£ - rfoC 1 - V l v C 2 - r,lC 3 = 0, 
V'CI + V 2 C 2 V + rfC 2 w + r?„ - X 2 rj 2 + (2tf - vl)C 2 - ^ - V 2 W C 3 = 0, 
V'CI + r, 2 C 3 v + r, 3 C 3 w + rf xx - X 3 T) 3 t + (2£ - rfJC 3 - ^C 1 - i? u C 2 = 0. 
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Because the functions C k (u, v, w) have the known structure defined in system (pQ) (otherwise 
the problem is very difficult even in the case of two-component RD system (7J IS]), this system 
of DEs can be solved in a straightforward way. In fact, solving subsystem (IBl)-( fT0l) . one obtains 



rj 1 = r l {t, x)u + q l {t)v + h}(t)w + x), 
rf = r 2 (t, x)v + q 2 (t, x)u + h 2 (t)w + p 2 (t, x), 
= r 3 (t, x)w + q s (t, x)u + h?(t)v + p 3 (t, x), 



(23) 



where £°, r h , q k , h k and p k are unknown functions at the moment. They can be found by 
substitution (1231) into ( fTTl) - (l2"2"l) and integration of the linear system obtained. The result will 
depend on the coefficients a^, dk and A& > 0, and can be formulated as theorem 1. 

Theorem 1 The DLV system (Op admits a non-trivial Lie algebra of symmetries if and only 
if one and the corresponding symmetry operator (s) have the forms listed in table 1. Any other 
DLV system admitting three- and higher-order Lie algebra is reduced to one of those from table 
1 by the local transformations: 



(24) 



u -> en exp(c 10 t)u + c 12 v + c 13 w, 
v C21 exp(c 20 t)v + c 22 u + c 23 w, 
w c 3 i exp(c 30 t)w + c 32 u + c 33 v, 
t — > c^t + C41, x — y cz>qx + C51, 

where Cij are some correctly- specified constants in each case (i = 1, . . . , 5, j = 0, . . . , 3/ 

Let us apply Definition 1 to find Q-conditional symmetries of the first type of the form ([3]) with 
(° ^ 0. In contrary to the standard Q-conditional (non-classical) symmetry, here we cannot 
simply set £° = 1 because some operators will be missed as it was shown in [16J. From the 
formal point of view, one needs to use three different manifolds in the criteria However, we 
take into account that the DLV system ([1]) has a symmetric structure and admits two discrete 
transformations of the form u v,v u,w w and u — >■ w,v — > v,w — > u. Thus, the 
algorithm is as follows. 

At the first step we apply Definition 1 in the case of the manifold {S\ = 0, S 2 = 0, S3 = 
0, Q{u) = 0} to find the operators ([3]) with £° 7^ 0. It turns out that the relevant calculations 
lead to a system of DEs with the very similar structure to (I8l)- (l22l . There are only two 
differences: no equations (TTTT) and the equations 

v'cl + v 2 c 2 v + v 3 c 2 w + (A X - x 2 )p 2 u + nl - X 2V 2 + (2g - vDc 2 - ^ - V 2 W C 3 = 0, 

rfCl + rfCl + r?Cl + (X, - A 3 )^ + rf xx - X^ + (2g - ^JC 3 - rf u C l - n\C 2 = 

are obtained instead of (12~T!)- (!2"2"|) . The simple analysis of these equations leads to the important 
conclusion: any DLV system (Q]j with Ai = A 2 = A3 admits only such Q-conditional symmetries 
of the first type, which coincide with Lie symmetries. Thus, we shall consider system (CQ) with 
at least two different lambda-s. 
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The second step of the algorithm consists in solving the obtained system of DEs. This step 
can be realized in a similar way as for the system of DEs (1B])- (122]) . however, all solutions leading 
to the Lie operators should be excluded (it means that the solutions obtained should not satisfy 
equations (llljl ). Moreover, all the Q-conditional symmetries obtained were reduced to the 
simplest forms using Proposition 1. Having this done, we arrived at exactly nine DLV systems 
with correctly-specified coefficients, which admit Q-conditional symmetries of the first type. 
The last step of the algorithm consists in applying Definition 1 to each of these systems using the 
manifold manifolds {Si = 0, S2 = 0, S3 = 0, Q{y) = 0} and {Si = 0, S2 = 0, S3 = 0, Q{w) = 0}. 
This allowed to construct all possible Q-conditional symmetry operators of the first type. The 
final result is formulated as follows. 

Theorem 2 The DLV system (T7]j is invariant under Q- conditional operator (s) of the first 
type (T3j] (with £° 7^ 0) if and only if one and the corresponding operator (s) have the forms 
listed in table 2. Any other DLV system admitting a Q-conditional operator of the first type 
is reduced to one of those from table 2 by a local transformation from the set (2~4\ ). Simulta- 
neously this Q-conditional operator is transformed to the corresponding operator listed in table 
2 (up to equivalent representations generated by adding a Lie symmetry operator of the form 
h 2 (t, x, u, v, w)d v + h 3 (t, x, u, v, w)d w ). 
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Table 1. Lie symmetry operators of the DLV system (OQ) 





Reaction terms 


Restrictions 


Additional Lie symmetries 


1 


u{b\u + c\v + d\w) 
v(b 2 u + c 2 v + d 2 w) 
w(b 3 u + c 3 v + d 3 w) 




D = 2i<9 4 + xd x — 2(ud u + vd v + 


2 


u(c\v + diw) 
v(a 2 + c 2 v + it;) 
w(a 3 + v + d 3 u>) 






3 


n(cif + d\w) 
v(c 2 v + u>) 
■u;(t> + d^w) 




u<9 u , D 


4 


u(a\ +bu + v) 
v(a 2 + u + cu ) 
+ cu) 


A2 = A3 = 1 


exp(-a 2 t)vd w , wd w 


5 


+ u ) 
w(u + cv ) 
+ cu) 


A2 = A3 = 1 


vd w , wd w , D 


6 


-u(ai + -u + f ) 
v(a 2 + u + v) 
w(u + v ) 


Ai = A2 = A3 = 1, 
a\a 2 (ai - a 2 ) 7^ 


exp(—ait)ud w , 
wd w , exp(-a 2 t)vd w , 
(a 2 (u + ai) + a\v)d w 


7 


it(a + u + 
w(u + 1>) 
w(u + t>) 


Ai = A2 = A3 = 1, 
a ^ 


exp(—at)ud w , 

wd w , vd w , 
(u + a + avt)d w 


8 


u(bu + 

+ cv) 
w(bu + ci>) 


Ai = A2 = A3 = 1, 
6^1, c^l 


D, wd w , ((b - l)u + (1 - c)v)d w 
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Table 2. Q-conditional symmetries of the first type of the DLV system (Op) 





Reaction terms 


Restrictions 


Q-conditional symmetry operators 


1 


u{a\ +bu + bv + dw) 
v(d2 + bu + bv + dw) 
w(a,3 + u + v + d$w) 


a\ ^ a 2 , 
(b-l) 2 + (d-d 3 ) 2 ^0 


^ + a 1- aM<% ftj 


2 


u{a\ + u + v + w) 
v(d2 + u + v + w) 

U)(d3 + U + V + w) 


(ai - a 2 ) 2 + (ai - a 3 ) 2 + 


Q 2 , i = l,...,6 


3 


u{d\ + u + v + w) 

V[d2 -\~ U -\- V -\- W ) 
w(ds + U + V + w) 


(A 2 - A 3 )ai - A 2 a 3 + A 3 a 2 = 0, 


Qf, i = l,..., 6, 

R, 1 Rwr>( a2 ~ a ^+\ii(ri ft \ 

Ot -+- pexp i x 2 -x 3 l j u \°v o w ) 


4 


u(d\ + u + v + w) 

v(d2 + U + V + w) 
w(dz + U + V + w) 


(A 2 - A 3 )ai - (Ai - A 3 )a 2 + 
(Ai - A 2 )a 3 = 0, 
(ai -a 2 ) 2 + «V0 


Qf, t = l,...,6 


5 


u(d\ +bu + v) 

v(d2 + U + Cv) 

w(bu + v) 


(6-l) 2 + (c-l) 2 /0 




6 


u(d\ + u + v) 

v(d2 + U + v) 

w(u + v) 




Qf, Qf, i = l,...,4 


7 


u(d\ + bu + cv) 

V(d2 + U + V) 

w(bu + i>) 


A 2 = A 3 = 1, 
ai(l — 6) = a 2 6(l — c) 


ft + ((1 - 6)u + (1 - c)v + a 2 (l - c))9 w 


8 


u(a + 6u + cv) 
v(a + u + v) 
w(bu + v) 


A 2 = A 3 = 1, 
6^1, c^l, 
6(2 - c) = 1 


ft + (i - c)d w + 

((1 - b)u + {1 - c)v)ip 4 (t)d w 


9 


u(di + u + v) 
v(d 2 + u + v) 
w(u + v) 


A 2 = A 3 = 1 


Ql i = l,...,5 
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In table 2, the following designations are introduced: 
Qi = Qi wrt h oc = 0, i — 1, ■ ■ ■ , 6; 

Qt = 9 t + ^ — ^-u(d u - d v ) + au{d v - d w ), Q\ = d t + — °^v(d v - d u ) + av{d u - d w ), 

Ai — A2 Ai — A2 

Q3 = <9t + ^ — - d w ) + ««(<% - d w ), Q\ = d t + — ^ — - <9 U ) + aw(d u - 

Ai — A3 Ai — A3 

Q5 = <k + v — T- V ( d v - + - d w ), Ql = d t + — ^-w(d w - d v ) + cra;(<9 u - d v ); 
A2 — A3 A2 — A3 

Ql = d t + atid x + exp (i^—^—al - a x )^- + — -aix)u<9 to ; 

Q!j> = a t + ai 3 x + exp ^ 2 a\ - a 2 )y + — - ^ 3 a 1 x^jvd w , 

Ql = d t + - d v ) + /3exp ^-^-(^ -^ U 

Ai — A 2 V A3(^A2 — AiJ / 

ft + ^f^(* - d u ) + /3ex P ( (A2 ~ A ^!~_ (A A 1 ~ A3)a2 t)^, 

12 3 v 1 2 / 

^ 6 a , a 2 Ai-aiA 2 /(A 3 - A 2 )oi - (A 3 - Ai)a 2 \ , , 
^4 = d t + t 7T r-rtu9„, + exp — — t)w{d u - d v ); 

A3\A2 — M) v As^Ai — A2J / 

Q\ = Q\ with A 3 = 1, Ql = Qt with A 2 = A 3 = 1, 

Ql = dt + -Y — ^r«(<9« - dv) + (<Pi(t)u + <p 2 (t)v + /3i)d w , 

Ax — 1 

Ql = d t + (tp 3 (t)u + ip 2 (t)v + /3x)d w , Ql = d t + ^ — ^-v(d v - d u ); 

Ai — 1 

where the functions y?j(t) (i = 1, . . . , 4): 



ftt + if a 2 = 0, = f ft*, if a 2 = 0, 

/3 2 exp(-a 2 t) + g, ifa 2 ^0; ^ U j if a 2 ^ 0; 



f . = hit + ^, ifai = 0, f£) f* + A ifa = 0, 

V ' 3U |/3 2 exp(-ait) + ^, ifai^O; ^ U \/?exp(-at) + I if a ^ 

while a and /3 (with and without subscripts 1 and 2 ) are arbitrary constants. 

Finally, we point out for the reader that (i) the inequalities listed in the second column of 
table 2 guarantee that the relevant operators from the third column are not equivalent to any 
Lie symmetry operators listed in table 1; (ii) if the given DLV system from table 2 contains 
as a particular case (up to local transformations of the form (|24l) ) a simpler system listed in 
another case of table 2 then, in order to get the complete list of symmetries, one should go the 
case with the simpler system . 
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Example. Let us apply the results obtained to a biologically motivated system of the form 
(CTJ. One notes that the DLV system in case 4 of table 2 is equivalent (the relevant substitution 
is u — y —bu, v — > —cv, w — >■ —dw) to the system 

X\u t = u xx + u(ai — bu — cv — dw), 

X 2 v t = v xx + v(a 2 — bu — cv — dw), (25) 
X 3 w t = w xx + w(a 3 — bu — cv — dw), 

where the coefficients at, b, c and d are the known positive constants. It is well-known that 
system ( 125]) is used for modeling competition between three species in population dynamics 
[2j H]. Substituting u — > —bu, v — > —cv, w — > —dw into the Q-conditional symmetry operator 
Q\ (see case 4 of table 2) one obtains 

Qi ->■ Q = d t + ^ — ^ m(<9 u - - <%) + a&«(- <9„ - ^ d w ) . (26) 
Ai — A2 c c a 7 

Applying the standard procedure for reducing the given PDE system to a ODE system via the 
known symmetry operator (T26]) . we easily find the ansatz 



' ai 


-02 




-A 2 


Ai- 


-A 2 


ai- 


-a 2 


Ai- 


-A 2( 



(27) 



. Ai— A 2 



where y?i(a;), V^^) and ^3(2;) are new unknown functions. Substituting ansatz (j27|) into (T25l) 
and taking into account the restriction (A2 — As)ai — (Ai — X^)a2 + (Ai — X 2 )a 3 = (see case 4 
of table 2), one obtains reduced system of ODEs 

¥i + Vi( Xl Z-x 2 ai -<P2-<p 3 ) =0, 

(p' 2 ' + ip 2 (a 2 - ¥2 - <fs) = 0, (28) 
ipl + ip 3 (a 3 -ip 2 - y? 3 ) = 0. 

Now exact solutions of the 3-D competition system ( !25l) can be easily derived by inserting 
solutions of this system into ansatz ( f2~Tl) . For example, the constant solution ip 2 = vq, tp 3 = 
02 — ^0 of the second and third equations of ( |28l) with a 2 = a 3 , produces the following exact 
solution of the 3-D competition system fT25|) with a± 7^ 0,2 = a 3 : 



u 



ift(»)««P(»). 

f + ^^S-^W^l^). (29) 



where y?i(a;) is a solution of the linear ODE <p'{ + A2 ^"1 A " 1 ' > <fi = 0. It should be stressed that this 



solution is not obtainable by any Lie symmetry because system (12511 is invariant only under the 
principal algebra ([6]), so that only traveling wave solutions can be constructed. 
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4 Conclusions 



In this communication Lie and Q-conditional symmetries for the 3-D diffusive Lotka-Volterra 
system ([I]) are studied. The main results are presented in theorems 1 and 2 giving the exhaustive 
list of the systems admitting Lie symmetry and Q-conditional symmetry of the first type, 
respectively. To the best of our knowledge, there is only paper [25] (see section 4.1) where 
the authors attempted to construct Q- conditional symmetry operators for a three- component 
PDE system, the classical Prandtl system, using so called the iterating non-classical method. 
However, all the results obtained therein are obtainable by the Lie method (this is clearly 
indicated in [25]). Here the Q-conditional symmetry operators listed in table 2 are not reducible 
to Lie operators from table 2. Thus, we have constructed non-Lie symmetries for a multi- 
component non-linear system of PDEs for the first time. Moreover, we have applied the Q- 
conditional symmetry operator ( 1261) to reduce the relevant Lotka-Volterra system and shown 
that even the simplest particular solution of the ODE system obtained leads to the exact 
solution (I27p . which cannot be derived by any Lie symmetry. 
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